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Abstract 



We consider models for the plateau transition in the integer quantum Hall 

effect. Starting from the network model, we construct a mapping to the Dirac 

Hamiltonian in two dimensions. In the general case, the Dirac Hamiltonian 

has randomness in the mass, the scalar potential and the vector potential. 

Separately, we show that the network model can also be associated with a 

nearest neighbour, tight-binding Hamiltonian. 
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I. INTRODUCTION 



Anderson localisation is central to understanding of the integer quantum Hall effect 
(IQHE) [|T| . In particular, the plateau transitions, between different quantised values for the 
Hall conductance, reflect delocalisation transitions in each Landau level. Scaling ideas |2| 
provide a framework for understanding these transitions, and are supported by the results 
of experiment [0] and of numerical simulation UJ. Progress towards an analytical theory of 
the critical point, however, remains limited. 

The simplest starting point for such a theory is to neglect electron-electron interactions 
and consider a single particle moving in a magnetic field with a disordered impurity poten- 
tial. In pioneering work, Pruisken and collaborators Q obtained from this a field-theoretic 
description, in terms of a u-model. More recently, in response to the difficulties of extracting 
quantitative results from the a-model, several alternative formulations have been explored: 
Read |J, Lee || and Zirnbauer have investigated spin chains; Lee and Wang |8j have con- 
sidered the replica limit of Hubbard chains; and Ludwig and collaborators || have discussed 
the Dirac equation. 

The correspondence between Dirac fermions in two space dimensions, and non-relativistic 
charged particles moving in a magnetic field, stems from the fact that time-reversal symmetry 
is broken both by a mass term in the two-dimensional Dirac equation 0,0, and by a 
magnetic field in the Schrodinger equation. Moreover, as emphasised by Ludwig et al, 
the Hall conductance of Dirac fermions, with fixed Fermi energy, has a jump of e 2 /h, if 
the fermion mass is tuned through zero. The critical behaviour at this transition depends 
on the symmetries of the Hamiltonian. The Dirac equation with only a random vector 
potential is particularly amenable to analysis |9,[lTl, since the zero-energy eigenstates are 
known explicitly ||12|| . Critical properties are controlled by a line of fixed points, and turn 
out to be different from those expected at plateau transitions in the IQHE. The line of 
fixed points, however, is unstable against additional randomness, either in the mass or in 
the scalar potential, and flow is conjectured M to be towards a generic quantum Hall fixed 
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point, describing the same critical behaviour as emerges from the usual Schrodinger equation. 

Confidence that Dirac fermions with suitable randomness do indeed have a critical point 
in the same universality class as the IQHE plateau transitions is clearly strengthened if there 
exists an explicit mapping from a microscopic model for the IQHE to the Dirac Hamiltonian. 
Fisher and Fradkin [TB[], and subsequent authors P,|H)[, have reached the Dirac equation 



starting from certain, rather specific, tight-binding models. An alternative to the tight- 



binding model, as a description of the IQHE, is the network model ||16|| , studied extensively 
by numerical simulation JlTj . Ludwig and collaborators || have asserted that Dirac fermions 
with the various possible kinds of randomness each represent particular forms of network 
model. These authors, however, did not set out a transformation from one model to the 
other. Separately, Lee found such a transformation, in the particular case of a network 
model without random phases, obtaining Dirac fermions with randomness only in the mass. 

The purpose of this paper is to describe a general mapping from the network model to 
the Dirac Hamiltonian in two dimensions, which, in the unrestricted case, has randomness in 
the mass, the scalar potential and the vector potential. Any approach to this problem must 
confront the fact that the network model is defined using the language of scattering theory, 
and therefore, at least in the first instance, contains information only about behaviour at one 
energy. The Dirac Hamiltonian, by contrast, obviously fixes properties of an entire spectrum 



of eigenstates. We begin from a unitary matrix defined |19| for the network model, which, 
heuristically, can be thought of as a time-evolution operator. We show, in a continuum 
limit, that it is the evolution operator for a Dirac Hamiltonian. In this respect, our route 
is rather different from that of Lee ||, who obtains a Hamiltonian by endowing the phases 
of the network model with an energy dependence. We also differ in taking the continuum 
limit isotropically, while Lee || does so anisotropically. 

Our mapping is described in section II. In section III we examine in detail how edge states 
of the network model are related to boundary states of Dirac fermions. This is important, 
since it is these states which are responsible for the quantised Hall conductance away from 
plateau transitions. 



Equivalence between the network model and the Dirac Hamiltonian necessarily requires 
a continuum limit. In section IV, we show that, independently of the continuum limit, 
one can associate with the network model a tight-binding Hamiltonian, which contains only 
nearest-neighbour hopping. 



II. MAPPING FROM THE NETWORK MODEL TO THE DIRAC 

HAMILTONIAN 



In this section we construct an explicit mapping from the network model JT6| to the 
Dirac Hamiltonian in two dimensions. First, we recall the physical basis for the network 
model and its definition. Consider non-relativistic, charged particles moving in a smoothly 
varying scalar potential in two dimensions, with a strong perpendicular magnetic field. 
The potential is smooth if its correlation length is much larger than the cyclotron radius, 
and the field is strong if the cyclotron energy is larger than the amplitude of potential 
fluctuations. Under these conditions, the kinetic energy of cyclotron motion about the 
guiding centre, and the potential energy associated with the position of the guiding centre, 
are both separately conserved. We focus on drift of guiding centres along equipotential lines. 
In the network model, portions of a given equipotential are represented by directed 'links', 
and the wavefuntion for the particle is caricatured by complex current amplitudes, Z, defined 
at points on each link. On traversing a link, a particle aquires an Aharonov-Bohm phase: if 
Zi and Zj are amplitudes at opposite ends of the link k (see Fig. la), Zj = e lr ^ k Zi. Tunneling 
between two disjoint portions of equipotential can occur where they are separated by less 
than a cyclotron radius, as happens near saddle-points in the potential. It is incorporated 
into the model at 'nodes', where two incoming and two outgoing links meet. The amplitudes 
on the four links that meet at a given node may be related by a transfer matrix or by a 
scattering matrix. In a suitable gauge, each of these 2x2 matrices is real and depends on 
a single parameter, which we denote by 9 (for the transfer matrix) and j3 (for the scattering 
matrix). The parameter determines the relative probabilities for a particle to turn to the 
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left or to the right on arriving at the node. It is a smooth function of the equipotential 



energy, measured relative to the potential at the saddle-point flIB| . Referring to Fig. (lb), 
one has 
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The two parameters are related by sin/5 = — tanh#. On varying the equipotential energy 
from far below that of the saddlepoint to far above, P increases from P = to P = 7r/2; 
tunneling is a maximum at P c = 7r/4. 

The network model as a whole is built by connecting these two elements - links and 
nodes - to form a lattice. The simplest choice is the square lattice, illustrated in Fig. 2. 
Randomness is introduced by choosing each link phase, <pk, independently from a probability 
distribution. The model represents particle motion at an energy determined by the value 
of the node parameters. If all nodes are identical, and if phases are uniformly distributed 
between and 2n, the system is critical at P — p c , and in the localised phase otherwise. 

We follow Klesse and Metzler |[19|| , and associate a unitary matrix with the model. 
Roughly speaking, this matrix is a time evolution operator. Let the unit of time be the 
interval required for a guiding centre to drift from the mid-point of one link, through a 
node, to the mid-point of the next link; ignore dispersion in this time interval, arising from 
variations in drift velocity or in lengths of links. Let Z(r\ L) be the amplitude for a particle 
to arrive at a point, r, after L time-steps, starting from an initial wavefunction Z(r';0). 
Then 



Z(t;L + 1) =J2TryZ(r';L), 



(3) 



and T is the required time evolution operator. Eigenf unctions of T with eigenvalue 1 are 
stationary states of the network model. 



In Eq |3], the element T r y is non-zero only if there is a one-step path on the lattice from 
r to r': that is, a path which follows the directions of the links and passes only one node. 
The values of these non-zero elements are given |L9| by a product of a phase factor from the 
link traversed, and a tunneling amplitude from the node, with sign conventions indicated in 
Fig. 3. 

To be definite, consider the system illustrated in Fig. 4. Plaquettes are labelled by the 
coordinates, (x, y), of their centres. With our choice of lattice constant and of orientation for 
the axes, (x, y) are a pair of integers, either both even or both odd. We denote the four links, 
i, making up a plaquette by i — 1, 2, 3 and 4, so that a point, r, on the network is specified 
by the combination (x,y,i). Initially, we take the tunneling parameter, (3, to be the same 
at every node, and the four phases, (pi, to be the same on every plaquette. In addition, it is 
convenient to measure the phases relative to their value when half a flux quantum threads 
each plaquette, by replacing 4 with 4 + ir. 

Because each plaquette has four links, the matrix T has a 4 x 4 block structure, and 
because the mid-points of links form a bi-partite lattice, each such block consists of two 
2x2 blocks. To exhibit this structure, we arrange the amplitudes Z(r; L) = Zi(x,y;L) in 
the order (Z + ,Z_), with Z + (x,y) = (Zi(x,y), Z 3 (x,y)) and Z-(x,y) = {Z 2 {x, y), Z±(x, y)), 
suppressing the time, L, for clarity. In this basis, the evolution operator is 
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Here, we have introduced the abbreviations C = cos (3 and S = sin /3, and the translation 
operators, t x ± and t±, defined by their action, t x ±Zi{x,y) = Zi{x ± l,y) and t y ±Zi{x,y) = 
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Zi(x,y ± 1). The first row of M, for example, expresses the fact, illustrated in Fig. 4, that 
Z x {x,y;L + \) = Se** l Z 2 {x - 1, y + 1; L) + Ce^Z^x, y; L). 

In order to decouple Z + from Z_, we consider the two-step evolution operator, 



W = T 2 
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We may then deal just with the upper-left block, U = MN, in the matrix W, and the 
component-pair Z + . Since, at this stage, we are treating a system without disorder, U is di- 
agonalised by Fourier transform. We write its eigenvectors, u, as u T (x,y) = [v,w)e % ^ qxX+qyy ^ 
and find 
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where 
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and 7*, a* are the corresponding complex conjugates. The eigenvalues of U are e^ x ~ v ^ with 
a phase, Xi given by 



cos% = sin2/3cos(g x - A x ) cos(q y - Ay), 



(10) 



in which A x = (<pi — <f>s)/2 and A y = (0 4 — fa)/ 2- Setting j3 = (3 C + m/2 = 7r/4 + m/2, and 
taking —it < x < 7r, the range of allowed values for x has gaps around x — and x — i 71 " for 
m 7^ 0: x satisfies — n + |m| < x < — l m l or l m l < X < ^ ~~ l m l- This dispersion relation is 
illustrated in Fig. 5. Stationary states of the network model are characterised by x — V = 0. 
To extract from the unitary evolution operator a Dirac Hamiltonian, H, we write U = 



-iH 



and work in the continuum limit, in which H is small. Thus we expand around 



(lx,Qy) = (0,0), taking m and the link phases, <fij, to be small. To leading order, Eq.(|T0"D 
gives the spectrum 



X 2 = m 2 + (q x - A x ) 2 + (q y - A y ) 



(11) 



for small x- At an operator level, when t± and t± act on smooth functions we make the 
replacements t± = 1 ± d x and t\ = 1 ± d y . Then 
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Finally, to bring the Hamiltonian into a conventional form, we make a rotation in the two 
component space, setting H = GHG^ 1 , with 
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H = (p x - A x )a x + (p y - A y )(jy + ma z + V , (14) 
-id xi and similarly for p y , and we use the Pauli matrix representation 
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Now consider a network model with randomness. If the link phases and tunneling pa- 
rameter vary smoothly in space, one obtains in the continuum limit the Dirac Hamiltonian, 
Eq. ([14]) , with randomness in the vector potential, scalar potential and mass. Specifically, 
fluctuations in the vector potential, A, arise from randomness in the individual link phases, 
fluctuations in the scalar potential, V, come from variations in the total Aharonov-Bohm 
phase associated with each plaquette, and fluctuations in the mass, m, are present if the tun- 
neling parameter is not constant everywhere. The time-independent states of the network 
model correspond to the zero-energy states of the Dirac Hamiltonian. 

This mapping can also be carried through for generalisations of the network model. In 



particular, the two-dimensional model in which each link carries N channels |T6j is equivalent 



to the U(N) Dirac Hamiltonians investigated by Fradkin [[14 
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III. EDGE STATES IN THE NETWORK MODEL AND THE DIRAC EQUATION 



The edge of a sample is, of course, set by a scalar confining potential in the usual 
description of the IQHE, based on the Schrodinger equation. Dirac fermions, by contrast, 
are confined by a spatially dependent mass, as discussed by Ludwig et al M. In particular, 
chiral, zero-energy states of Dirac fermions are associated with contours of zero mass ||. 
We discuss in this section how such edge states emerge in our mapping from the network 
model to the Dirac Hamiltonian. 

Consider a network model defined on a strip of finite width, as in Fig. For energies in 
the lower half of the Landau level, corresponding to values of the node parameter (3 < /3 C , 
all states are localised, while for energies in the upper half of the Landau level, for which 



(3 > p c , a pair of extended states appears, one at each edge of the strip [|16[ . 

The evolution operator, U, acting on the two-component wavefunction, Z + (L), intro- 
duced for the bulk of the system in the preceeding section, is supplemented at the edge by 
the following boundary conditions (see Fig. |^): at x — 0, the component Z 3 obeys the same 
equation as in the bulk, 

Z 3 (0,y;L + 2) = [UZ + (L)} 3 (0,y), (16) 

while the component Z\ satisfies 

Z 1 {0,y;L+2) = -e t[M °> y)+M °' v)] Z 3 (0,y, L). (17) 

We wish to check under what conditions the evolution operator, U, has an eigenvector, u 
representing an edge state. We simplify the discussion by considering a semi-infinite system 
without disorder, setting fa = for all i. Without disorder, the spatial dependence, for 
x > 1, of such an eigenvector is u T (x,y) = (v,w)e^ qy ~ Xx \ where Re[X] > 0; for x = 0, 
u T (0,y) = (v',w')e iqy , where v,v',w and w' are constants. This ansatz in the equation 

Uu = e ix u } (18) 

taking U from Eqs. (0), ([IB]) and (|17|) , yields 



tan j3 



(19) 



confirming that an edge sate exists (i?e[A] > 0) only for /3 > /3 C = 7r/4. 

Similar results also follow if one considers directly the continuum limit. Let the eigenfunc- 
tions of H be \j>(x, y), so that those of H are ^(x, y), with \P = G^/. Writing \I/ T = (vE^, ^> 3 ), 



(20) 



the boundary condition at x — 0, Eq. (|17D, is to leading order 



* 3 (0,jO = itti(0,y). 



Note that this boundary condition enforces a chiral edge current: the current density, with 
components j a = ^Vq,^, is j = (0, |^i| 2 ), and necessarily in the positive y-direction. Im- 



posing this boundary condition, the Dirac Hamiltonian, H, of Eq. (|14| ) with V = and 
A = has an eigenfunction of energy E 

(A 



V(x,y) = e iEy e- mx 



(21) 



provided that m > 0. This same eigenstate appears from considering an infinite system 
with position-dependent mass, m(x,y), following Ludwig et ah setting m(x,y) = m for 
x > 0, and m(x, y) = m for x < 0, the boundary condition, Eq. (pop, emerges in the limit 
mo — > — oo. 



IV. MAPPING FROM THE NETWORK MODEL TO A TIGHT-BINDING 

MODEL 

It is also possible, without taking a continuum limit, to associate a nearest-neighbour 
tight-binding Hamiltonian with the network model. The sites of the tight-binding model, 
each carrying one basis state, correspond to the links of the network model. In terms of the 
one-step evolution matrix, T, the tight binding Hamiltonian, 7i, is simply 

U = (T f - T)/i. (22) 
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We indicate schematically in Fig. which matrix elements of TC are non-zero, and give their 
values in terms of link phases and the tunneling parameter. 

This Hamiltonian has two important features. First, it is natural to introduce a unit cell 
containing the four sites arising from one plaquette of the network model. The amplitudes 
of nearest neighbour hopping within and between unit cells have moduli cos (3 and sin/3, 
respectively, where (3 is the tunneling parameter: they are different, except at the critical 
point, (3 = (3 C = 7r/4. Second, the phases of the hopping matrix elements are correlated in 
the way indicated in Fig. 0. It follows from known behaviour of the network model that these 
correlations have unusual consequences for the tight-binding model. Consider a system in 
which all link phases are independently and uniformly distributed. It is straightforward to 
see that for this system the eigenvalues of T in the complex plane are distributed uniformly 
on the unit circle. As a result, the density, p(E), of eigenvalues, E, of TC can be given exactly: 
p(E) = 1/471- [l-(£/2) 2 ] 1/2 for (E/2) 2 < 1 and p(E) = 0for(£/2) 2 > 1. In addition, one can 
see that the eigenvectors of T have the same statistical properties throughout the spectrum, 
a feature inherited by TC. Hence the localisation length of eigenstates of TC is independent 
of their energy, E. If all nodes have the same parameter value, (3, then as (3 — > (3 C , the 
localisation length diverges uniformly across the spectrum: this Hamiltonian never has a 
mobility edge function of energy. 

We note that TC is similar in structure to, but different in detail from the tight-binding 
model used as a depature point by Ludwig and collaborators ||. The latter includes not 
only nearest-neighbour, but also next-nearest neighbour hopping: compare Fig. |7| with Fig. 



1 of Ref [91. Our model also differs from that of Fisher and Fradkin [13 



V. SUMMARY 

We have set out in detail a mapping, from the network model for plateau transitions 
in the IQHE, to Dirac fermions in two space dimensions. The mapping makes crucial use 
of a unitary operator defined for the network model |19||, which is essentially the time- 
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evolution operator. The two-component structure of Dirac spinors in two space dimensions 
arises rather naturally from the network model, defined on a square lattice: the fact that 
each plaquette has four sides suggests a four-component wavefunction, which separates into 
two independent pairs because of the existence of two sublattices. This structure is not 
dependent on the continuum limit, and is also shared by a nearest-neighbour tight-binding 
Hamiltonian, derived directly from the evolution operator. 
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FIGURES 



FIG. 1. Components of the network model: links (a), and nodes (b). 
FIG. 2. The network model on a square lattice. 

FIG. 3. Amplitudes associated with possible scattering paths at nodes. 

FIG. 4. The network model, showing: (i) our coordinate system for plaquttes; (ii) labelling of 
the four links that make up a plaquette; and (hi), with dashed lines, the paths that contribute to 
Zi(x, y;L + l) for i = 1,3. 

FIG. 5. The dispersion of x ( vertical axis ) plotted as a function of (p — A) x and (p — A) y 
(horizontal plane) in units of n. The width of the gap between the two bands is determined by the 
mass which is 1.2 here. 

FIG. 6. The network model defined on a strip of width M. Dashed arrows indicate the propa- 
gation direction of edge states. The dotted arrow represents the boundary condition, Eq 17. 

FIG. 7. Schematic illustration of the tight-binding Hamiltonian. (a) The only non-zero ma- 
trix elements are those linking nearest-neighbour sites, within plaquttes (full lines) and between 
plaquettes (dashed lines), (b) Their values are: %e l ^ x - y ^ cos/3 for the bonds marked I = 1,2,3,4; 
i e i<f>i(x.y) s [ n p^ f or bonds marked I = l',2'; and —i^i^-v) sin/3, for the bonds marked I = 3', 4', 
hopping in all cases being in the direction given by the arrows. 
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